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^ ' Abstract 

o ■ 

O^' The studies of influence of spin on a photon's motion in a Schwartzschild spacetime 

I , is continued. In the previous paper the first order correction to the geodesic 

motion is found for the first half of the photon world line. The system of equations 
for the first order correction to the geodesic motion is reduced to a non-uniform linear 
^ . ordinary differential equation. The equation obtained is solved by the standard method 

H ' of integration of the Green function. 

1 Introduction 

Existence of spin-gravitational interaction was proved by A. Papapetrou in his analysis of 
motion of deformable body [H |2] . More recently we have derived Papapetrou equations as 
a reduction of equations of motion of (tangent) rigid body [3]. However, analysis based on 
classical mechanics does not admit passage to the limit of zero mass. The problem of motion 
of a massless spinning particle is especially interesting due to importance of electromagnetic 
waves in astronomial observations. Importance of the problem of photon world line with 
account of the spin-gravitational interaction was under discussion for decades jl]-|Z!, however 
it could not be formulated properly until equations of motion for massless spinning particle 
were derived. 

Passage to classical mechanical limit from a non-scalar field equation is an alternative 
approach tried by numerous authors [H]-[ID1- Straightforward derivation of photon's equa- 
tions of motion from Lagrangian of electromagnetic field was made in our works fT|^]. As 
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the result we have obtained Papapetrou equation in the form 



uj-— = 



(1) 



where u is frequency and a"^ - spin of photon which hes on tangent vector subspace of 
polarization vectors fia, a, f3, ... = 2,3. The spin remains parallel to itself in the polarization 
subspace along the world line. It is seen from this equation that the effect of spin-gravitational 
interaction grows linearly with the wavelength, so, under large enough wavelengths it may 
well be observable. 

The next step is to find out photon world lines in a given model of space-time. The 
simplest way is to construct world lines via computing the first order correction to isotropic 
geodesic. These computation in case of Schwartzschild space-time was completed in our 
recent work ^2] where, however, the desired result was not reached because the method used 
is valid only on a half of the world line which starts from infinitely distant source and ends 
up at the minimal value of the radial coordinate r (periastr). To construct the whole world 
line one needs to match two such halves built in different coordinate systems that is quite 
difficult. The goal of the present work is to work out a method of building the whole of the 
world line of a massless spin 1 particle in Schwartzschild space-time. 

In the framework of the accepted approximation scheme we chose an isotropic geodesic 
(hereafter: reference geodesic) which is a world line of a spinless particle, then include spin 
into the right-hand side of the equation (jT)) putting it parallel to itself along the geodesic. 
It is convenient to introduce a geodetic flow in the equatorial plane (^^ = 7r/2 in spherical 
coordinates) as a vector field of velocities and its small variation. By construction the vector 
field satisfies geodesic equation, so it remains to find out its small variation such that the 
total vector satisfies the equation (^. This yields a linear equation of second order for the 
variation which can be solved on a geodesic of the flow. Solution of this equation is the first 
order correction to be found. The only independent variable in the equation is parameter s 
on the reference geodesic. 

The metric of the Schwarzschild space-time can be specified by a field of orthonormal 
frame. A sample of pair of dual to each other covector and vector standard (with no isotropic 
elements) orthonormal frames is 




(2) 



where orthonormality of the frames means 

< e' >= f]''", < ea, e, >= r]ab, e\e,) = 
rjab the standard Minkowski metric. The corresponding connection 1-form is 
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The 2-form of curvature referred to the frame is: 



1^2° = -^e^Ae^, f^s' = "-^e^AO^ (4) 

2 Isotropic geodesic flow, canonical parameter and the 
CO- moving frame 

We construct a congruence of isotropic geodesic hues lying wholly on 6* = 7r/2 hypersurface 
(an equatorial plane) by solving the Hamilton- Jacobi equation for isotropic geodesies on it 

(d* , d^) =0 (5) 

where the function to be found is of the form 

^ = Et-L^ + R(r) (6) 

and contains two constant parameters E and L. Accordingly to Hamilton- Jacobi theorem 
the integrated equations of geodesies are 

— — = const, —— = const. 
oE oL 

These equations specify the variables t and if as functions of the variable r: 

t = -^, ^ = ^0 + ^, 0^./2 (7) 
where the function R is found by integrating the Hamilton- Jacobi equation (0): 



pr ^/E'^r'^ — L^fl — r„/r) 

Rir) = -e{t) / ^ -JY^dr, (8) 

■Jro r[l — rg/r) 



where e{t) is the well-known step function equal +1 for positive t and -1 for negative t. 

Each geodesic of the congruence starts at r = oo under t = — oo and reaches the 
periastr r = ro at the moment t = 0. After this moment the cordinate r grows and reaches 
infinity under t —> -(-oo. We choose initial conditions such a way that (p\^^_^ = 0, so that 



^0 



OR 



.It must be noted that the isotropic geodesies of the congruence he wholly 

on surfaces of level of the function ^ or, in other words, \l/ is constant on each of them. 
Therefore this function cannot be used as a parameter on the geodesies. Instead, a canonic 
parameter on the geodesies will be introduced below. 

In order to reduce the number of non-zero components of spin of the particle it is 
convenient to introduce a co- moving frame on the geodesies like that done in work 
Since the geodesies are isotropic the co-moving frame contains isotropic elements n±. One of 



3 



them, n_ is orthogonal to surfaces \1/ 
tangent to them everywhere): 



const, (thus, due to properties of isotropic vectors. 



n_ = (dx\d^) d. 



and another isotropic vector n_|_ is chosen only due to the normalization condition (n_ 

= 1. Two other vectors which constitute polarization subspace n^, a = 2, 3 are normalized 

to unit, space-like, orthogonal to each other and to the isotropic vectors: 

(1 - Tg/r)-^ dt + Dr-^ - R'E-\l - Tg/r) dr, 

(1 - r,/r)/2 [(1 - r,/r)-i dt - Dr'^ + R' E-\l - Vg/r) dr] , 

-1 ^ (9) 

r '-Off, ^ ^ 



where 



{r^ - D^{l-rg/r)yh 



E! = e{t) 



'd - 



dR' 
dL 



(1 - rg/r) dr^ 



J^2^2 _ 1,2^1 _ Tg/r) 



r(l — rg/r) 

and D = L/E is the impact parameter. The vector n_ is obtained from the 1-form of 
momentum o?^, by index lifting therefore, it satisfies the geodesic equation: 

o = {n^ ■ V)n- = 

and hence n_ is the canonical velocity on the geodesic [15J. Thus, the canonical parameter 
can be found by equating n_ to i = dx^dsdi that gives: 

- = -R/Eil-rg/r), ^ = 

It is seen that sign of dr/ds changes when passing through the hypersurface t = which 
separates the isotropic geodesic into two halves. Under t < coordinate r decreases along 
the geodesic until at t = it reaches its minimal value tq and afterwards grows up to infinity 
under t > 0. Therefore r can be used as a parameter only on a half of the geodesic where 
the coordinate r runs the range r G (rg, oo). 

In order to compute the first order correction to the whole of the geodesic we need a 
monotonic parameter s on it. It is convenient to define it such that at the moment of time 
t = when the coordinate r takes its minimal value (r = Tq) the value of s is zero, for 
example: 

r dr 



e{t) 



(10) 



1^2 _ 1)2(^1 _ Tg/r) 

Now we can find the relation between two parameterizations on the geodesic for each of its 

e{s)rdr 



halves: 



ds 
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^r^ - D^{l~rg/r) 

Since hereafter we shall use the values of the function R on the geodesic it is convenient to 
introduce an alternating function R{s) with its derivative w.r.t L : 



dRjs) 
dL 



sis) 



dR{r) 



dL 



dR 



dL 



Ddr 



ro ^ 



D2(l - rg/r) 



(12) 
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Besides, we need expressions of vectors fip = n'^ea via vectors Ca of standard frame: 



n_ 
n+ 

^2 
^3 



(1 - rg/r)-i/2 eo + Dr^i eg - R'E-\l - rgjrfl'^ Ci, 



62, 



(13) 



3 First order approximation of the Papapetrou equa- 
tion 



Assuming that interaction of spin with gravitation is weak and deviation of the ray from the 
reference geodesic: bx"" jr <^ 1 is small we shall derive the equation for the deviation from 
the Papapetrou equation and solve it. It must be noted that analysis made in the work 
contains one more inaccuracy: when considering deviation from a geodesic one should take 
into account the difference between values of connection on the geodesic and on the deviation 
from it. If it is done properly instead of the l.h.s. of the geodesic equation n_ o n_ one has 
that of the Jacobi equation 

|14j . Therefore the desired equation for the small deviation 6x^ extracted form the Papa- 
petrou equation ((T)) is 



ds^ ^ 



-XR^^,ntn^,nia'', 



(14) 



where A = cj^^ = const is the (asymptotical) wavelength. 

It was shown in the work JH] that the world line of a photon lie initially in the equatorial 
plane 9 = 7i/2 the vector of deviation has in the frame ()13|) the only non-zero component 
Therefore it is convenient to have explicit form of the following contractions of the Riemann 
tensor: 



P 2 

^bc2 



n n 



3rgD^ 
2r5 ' 



(15) 
(16) 



The latter is needed because in the co-moving frame (|13p the spin of photon has only one 
non-zero component o"^^. Taking into account that in the co- moving frame all components of 
the connection 1-form are zero we substitute partial derivative of the deviation vector for its 
covariant derivative and, referring to (jlj) and (fT3|) we reduce the Papapetrou equation for the 
deviation vector in Schwarzschild background to the following ordinary differential equation 



d^y , ^TgD^ 

+ ^ . y 



ds^ 



2r5 



3XrgD 
2r5 
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where the function y{s) to be found stands for Sx"^. Accordingly to fact that under s — > — oo 
the ray coincides with reference geodesic sought solution of the eq. must obey asymptotical 
conditions as follows: 

y(-oo) = 0, ^(-oo)=0. (17) 

It is convenient to rewrite the approximated Papapetrou equation in the standard 
denotions 



with 

Q{s) 



^^,+Q{s)yis) = F{s), (18) 



2r^ 



DJr^ - D^n -rJr) 
F{s) = -C^ ^ (19) 

C = 3Xrga^y2. 

This is a well-known linear non-uniform ordinary differential equation of second order which 
can be solved by standard methods. To do it we need first the solution of the uniform equation 
which corresponds to commonplace deviation of geodesies, thus, to the Jacobi equation in 
which spin-gravitational is 'switched off' by putting C = 0. Then as the uniform equation 
is solved we can use the solution for constructing solution of the non-uniform equation and 
obtain deviation caused by the spin-gravitational interaction. 



4 Integration of the Jacobi equation 

In this section we make straightforward calculation of the deviation vector which in the field 
of local frames p3|) has only one non-zero component 5x'^. For this end we consider a geodesic 
flow with small 2-component found from the corresponding solution of the Hamilton- Jacobi 
equation. Unlike the solutions considered above this solution depends on all four coosdinates 

^ = Et- + R{r) + e{e) (20) 

and contains three arbitrary constants E, L and M. The function R[r) is the same as before 
and the new function Q{6) has the form 



0(9) = ±f - My sin^ e. 



If we put L = M theis function vanishes and the corresponding geodetic flow reduces to 
the flow considered above. So, if we put (L^ — M^)/L^ = e"^ an infinitesimal parameter we 
obtain a family of geodetic flows which are infinitesimally close to the geodetic flow in the 
equatorial plane. Applying the Hamilton- Jacobi theorem to the solution (j^UJ gives: 

- -V? =F / d9 = const, (21) 



Ld9 OR 

— = i / . : H = const 

dL J JL^-Mysin^e 9L 
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and all the rest equations remain unchanged. 

The second line specifies the component 59 of the deviation 

y/L^ - M2 /dR \ . fdR \ 

66 = sm — h const = esm \ — h const 

L \dL J \dL J 

and the first line contains the same integral, therefore, neglecting the second poweer of the 
parameter e we obtain the same equation as under e = 0: 

..,r dR 
—if + ML — — = const, 

oL 

consequently, the component 5(p of the deviation is zero. Thus, the deviation vector has only 
one non-zero component which is 59: 

5r = 5t = 5ip = 0, 59 = e sin + (j , 

where ( = const and e is an abitrary infinitesimal parameter. Parameter ( specifies deviation 
from the geodesic at t = — oo which can also be put non-zero. In our denotions y = 5x'^ 
solution of the Jacobi equation just obtained is: 

y = ersm{RL + (22) 
where |i?L(ri)| stands for the elliptical integral (fT^ 



5 Integration of the non-uniform equation 

In this section we solve the non-uniform equation (fTH|) by standard method of integrating the 
Green function. By definition the Green function of this ordinary differential equation (fTHjl 
is a piecewise-smooth function composed of two solutions of the uniform equation. Since 
the differential equation is of second order the relevant Green function G{s\si) should be 
continuous in the point s = Si wile its first derivative makes a unit jump [T^. This way we 
construct the desired Green function by matching two solutions of the uniform equation as 
follows: 

r( \ \ - Hl ^ s <si, . . 

^^'\'^) - D \ rsm{eis)\RLir)\-eis,)\RLin)\} if s, < s; ^^^^ 

Here the factor ri/D normalizes the jump of the derivative, ri = r(si), r = r{s) are specified 
by the function inverse to the function s{r) defined in the equation (fTUI) and the variable s 
enters only the step function e{s) which specifies only the sign depending on the half of the 
world fine. Besides, solution under s < si is trivial that provides validity of our asymptotical 
conditions (fTTjl at s ^ — oo. 

Now, following the standard procedure we can find out solution of the nonuniform 
Jacobi equation (fT^ as follows: 

oo 

y{s) = f G{s\s^)F{s,)ds, = 

(24) 

sm{e{s)\RLir)\-e{si)\RL{n)\} 



— oo 
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Our goal is to describe behavior of the solution obtained under s > corresponding to 
deviation the ray off the reference geodesic after the periastr. It is convenient to represent it 
by the angle of deviation: 

, = ^ = -Cr {ejs) \R,ir)\ - ejs,) \R,in)\ } _ ^^^^ _ ^^^^^^ ^ 

T J —oo 



u s 



where we divide the domain of integration into parts of definite signs of the step function of 



'rl-D^{l-rg/ri) '''' ^Jrf - D^{1 - Vg/n) 

Here the square roots in the denominators appear when passing from the varibale si to ri 
fjlll At the same time we substitute the limit = oo of ri for the limit — oo of Si. After 
cancelling the square roots the expression simplifies and takes the form 



ro Ti •' '"0 



Introducing the notion of asymptotical angle S6oo = 66{+oo) which describes the final angular 
divergence between the ray and its reference geodesic we obtain expression for the deviation 
on s > branche of the ray as follows: 



r 





/■oo 


|i?i(oo)| 


/ COS 




Jro 



Evidently, the last term vanishes under s — >■ oo. Now let's find expression for 66 oo- After 
some algebraic operations we obtain: 

/•nci 

[Ur^. (25) 

It is seen that the integral \RL{ri) \ cannot be expressed in terms of convenient functions. 
However, since we are seeking for the first order on r^/ro -C 1 and the constant factor C = 
3rga'^^/2E already contains the factor Vg, it suffices to take the integral in zeroth order 
approximation on Vg/rQ. In other words, to obtain the desired result we can put Vg = 
in |i?i(ri)|. This is possible because though the integrand has square root singularity at 
ri = To all the integrals in (j25p are convergent. Moreover, the integrals as functions of Vg are 
continuous and differentiable at Vg = 0. 

Now, simplifying the integral (fT^ we obtain 

\RL{r)l.^^Q = arccos(ro/r). 

In this approximation we have \Rl{c>o)\ = n/2. Substituting these expressions into 
simplifies the integration and finally we have: 

dri f°° ro dri r„a'^^ Tqcy'^^ 



-2C / cos[arccos(ro/ri)J— - = — 2C / — ^ 



Erl ED^ 
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where we take into account that ro|,,. = D. Substituting this estimation to formula for the 
deviation we obtain under s > following asymptotical presentation: 

y{s) = r6e^ + ^6e^ + OiCr-^). 

This expression is valid for the second half of the world line and represents photon world 
line under large r. Since the deviation on the second half grows linearly even if there is no 
interaction, just because the initial plane of the ray is lost during the first half, this linear 
growth yields about one order grater deviation at infinity than that on the periastr obtained 
in the work 
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